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Abstract We report our observations on 4 AGNs viz, 
Mkn421, Mkn501, 1ES1426+428 and ON231 belonging 
to a sub-class called blazars. The observations were car- 
ried out using the Pachmarhi Array of Cherenkov Tele- 
scopes and span about 6 years period from 2000 to 2005. 
We discuss our methods of analysis adopted to extract 
the gamma ray signal from cosmic ray background. We 
present our results on the emission of TeV gamma rays 
from these objects. Also, we report on the status of the 
new high altitude experiment, a 7 telescope array at 
Hanle in the Himalayas at an altitude of about 4200 m 
above mean sea level for detecting celestial gamma rays. 

Keywords Mkn 421 • Blazar • 7-ray astronomy 



1 Introduction 

AGNs have dominated extragalactic 7-ray astronomy by 
virtue of their great luminosities. The general under- 
standing is that a supermassive black hole of mass 10 6 — 
10 9 M at the center of AGN accretes mass from the sur- 
rounding medium, forming an accretion disk and two jets 
emanating perpendicular to the plane of the accretion 
disk pp. These jets channel a plasma flowing out with rel- 
ativistic speed and any radiation produced inside them 
is greatly modified by Doppler effect. AGNs with jets 
directed towards us are called blazars. Blazars are char- 
acterized by two distinct parts in their Spectral Energy 
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Distributions (SEDs). First part in SED rises smoothly 
from radio wavelengths upto a broad peak spanning the 
range from optical to X-ray wavelengths and is due to 
relativistic electrons radiating via synchrotron process. 
Second part is probably due to inverse compton scatter- 
ing of synchrotron photons by the same electrons and 
is characterized by a peak in SEDs in hard X-ray to 7- 
ray band [5] • One of the characteristic features of these 
blazars is their time variability on scales ranging from 
minutes to years. Mkn 421, Mkn 501, 1ES1426+428, 
ON231 are the four blazars observed using Pachmarhi 
Array of Cherenkov Telescopes (PACT) during 2000 to 
2005. In the following sections we present our observa- 
tions, analysis procedure and the results obtained. Also, 
we have attempted a comparative study of two states of 
Mkn 421 namely, flaring and quiscent states, using radio, 
optical, X-ray and 7-ray data. Our future plans for the 
observations of 7-ray sources using a high altitude array 
are also presented. 



2 Pachmarhi Array of Cherenkov Telescopes 

Pachmarhi Array of Cherenkov Telescopes (PACT) is 
located in Central India (latitude 22° 28' N, longitude 
78° 25' E, altitude 1075 to). We use wavefront sam- 
pling technique to detect TeV 7-rays from astronomi- 
cal sources. There are 24 telescopes spread over an area 
of 80to x 100m. Figure [1] shows the schematic layout 
of PACT. Each telescope has 7 para-axially mounted 
parabolic mirrors of diameter 0.9m with a PMT (EMI 
9807B) at the focus of each mirror as shown in figure 
[2] Entire array is sub-divided into 4 sectors with 6 tele- 
scopes in each. Each sector has its own data acquisition 
system (DAQ) where data on real time, relative arrival 
time of PMT pulses (using TDCs) and photon density 
(using ADCs) of six peripheral mirrors in a telescope are 
recorded. Apart from this, there is also a Master DAQ 
at the center of the array for recording information of an 
event relevant to entire array. PMT pulses of 7 mirrors in 
a telescope are linearly added to form a telescope pulse 



2 



D. Bose et al. 



Table 1 Energy thresholds and collection areas for PACT 
at different incident angles 



Inclination 
Angle 


Energy 

Threshold 

(TeV) 


Collection 

Area 

(m 2 ) 


Vertical (0°) 


0.75 


1.4xl0 5 


15° 


0.90 


1.5xl0 5 


30° 


1.20 


1.8x10 s 


45° 


2.20 


2.7xl0 5 




Fig. 1 A Layout of PACT. The big circles represent the 
telescopes. Seven smaller circles inside a big circle represent 7 
mirrors in a telescope. Rectangular boxes represent the data 
acquisition centers. 

for trigger generation. Data recording is initiated when 
a coincidence of 4 out of 6 telescope pulses generates an 
event trigger for a sector. The typical trigger rate was 
about 2-3 Hz per sector. The orientations of telescopes 
are controlled remotely and monitored throughout the 
observations [3J. Details of the setup can be found else- 
where [HE]. To estimate energy threshold, collection area 
etc, we have carried out Monte Carlo simulations of ex- 
tensive air showers using CORSfKA package developed 
by KASKADE group 0. Energy threshold of PACT is 
estimated to be 750 GeV for vertically incident show- 
ers initiated by 7-raysQ and the corresponding collection 
area is 1.38xl0 5 m 2 . The flux level for detecting 7-rays 
from a source at 5cr sensitivity in 50 hours of obser- 
vation using PACT is estimated to be 2.93xl0 -11 ph 
cm s^ 1 assuming no cosmic ray rejection. In terms 
of Crab nebula flux this corresponds to about 0.9 Crab 
units above energy threshold of 750 GeV. For inclined 
showers, simulations were carried out using IACT op- 
tion in CORSIKA. The energy threshold and collection 
area increases with incident angle and these parameters 
are summarized in table [T] 

1 Assumed energy spectrum F(E)=kE _2,4 dE 




Fig. 2 Picture of one telescope with seven parabolic mirrors 
and PMTs 

Table 2 Observation log for blazars 



Year Duration of Observations (minutes) 





Mkn 421 


Mkn 501 


ON 231 


1ES1426+428 




z=0.030 


z=0.034 


z=0.102 


z=0.129 


2000 


3510. 


710. 






2001 


1960. 








2002 


1860. 


510. 




1520. 


2003 


1770. 


840. 


510. 


570. 


2004 


2270. 


780. 


550. 


870. 


2005 


930. 






960. 



3 Observations and Data Analysis 

Observations were carried out using PACT on clear moon 
less nights. Observations on source (source runs) were 
usually taken by pointing all telescopes to the source di- 
rection. The typical run span was about 1-3 hours. Back- 
ground runs were taken either immediately before or af- 
ter (sometimes both before and after) the source run by 
aligning all telescopes to a dark region (a region with the 
same declination as that of the source but with different 
RA). Background region is chosen in such a way that 
it covers same zenith angle range as that of the source. 
The observation log is given in table [H For each source 
an equal amount of data were collected on corresponding 
background runs. 

PACT data were analysed in the following way. Ce- 
lestial 7-rays are not affected by the interstellar magnetic 
field, therefore they retain their directionality. Whereas 
cosmic rays, being charged particles, are scattered by the 
interstellar magnetic field, as a result they are isotropic. 
Thus a source emitting 7-rays will be reflected as an ex- 
cess of events from the source direction compared to off- 
source direction. A number of preliminary checks were 
carried out on the data before doing actual analysis. A 
cut is imposed on the number of telescopes with valid 
TDC data to be > 8. Care is taken to see that the dis- 
tribution of telescopes with valid TDC data are simi- 
lar in both source and background data sets. The ar- 
rival direction of each shower is determined by recon- 
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Background normalized in 2.5 — 6.5 deg. 
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Fig. 3 Typical space-angle distribution of a source run 
(solid), background run (dotted) and the difference between 
the two (dashed). 



structing shower front using the relative arrival times 
of Cherenkov photons at various telescopes (or PMTs). 
Cherenkov photon front is then fitted with a plane, nor- 
mal to this plane gives the direction of the shower axis. 
Then, for each shower or event, the space angle is esti- 
mated as an angle between the direction of shower axis 
and the source direction. Thus space angles arc obtained 
for all events in source as well as background runs. Space 
angle distributions of all source runs are compared with 
the corresponding distributions of background runs over 
the same zenith angle coverage. Figure [3] shows the space 
angle distributions of events from source and background 
runs taken in a night. Space angle distribution of back- 
ground events is normalised to the source distribution by 
comparing the shape of the distributions in 2.5° to 6.5° 
window since we do not expect any 7-ray event in this 
region 5 . This normalisation is necessary since there are 
variations in the sky conditions at different times of the 
same night. Differences between the number of source 
and background events is calculated for each bin as (SV 
cBi) where c is a constant. We define, 

6.5 



J2 (Si - cB z ) 2 



i=2.b 

and normalisation constant c is chosen such that x 2 is 
minimum. The difference between the source and nor- 
malised background events in 0° to 2.5° region is then 
used as 7-ray signal. Thus time averaged 7-ray signal is 
obtained for each night observations. 



3.1 Analysis of Multiwavelength Data of Mkn 421 

We have attempted a comparative study of high and low 
states of Mkn 421. We chose a flare in March/April 2001 
as representative of high state as this flare is one of the 
strongest flares and February/March 2003 as a repre- 
sentative of low state. We have analysed X-ray archival 
data from Proportional Counter Array (PCA) on board 
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Fig. 4 2-9 keV light curve obtained from PCA data, upper 
panel for 2001 and lower panel for 2003. 



RXTE obtained during these two periods. The first data 
set corresponds to the period 19th March to 1st April 
2001 when Mkn421 was in flaring state. The second set 
was collected during 26th February to 5th March 2003, 
Mkn421 was in quiescent state at this time. PCA data 
were extracted from archival data sets of RXTE satellite 
from observation ids 60145 and 80172 respectively. 

The PCA consists of five identical xenon filled pro- 
portional counter units (PCUs) covering an energy range 
of 2-60 keV. During these observations only PCU0 and 
PCU2 were used. We have analyzed Standard 2 PCA 
data which has a time resolution of 16s with energy infor- 
mation in 128 channels. Even though observations were 
carried out with PCU and PCU 2, we have used only 
PCU 2 data because PCU had lost its front veto layer 
at the beginning of year 2000. So the data from PCU0 
are more prone to contamination by events caused by 
low-energy electrons entering the detector. Data reduc- 
tion is done with FTOOLS (version 5.3.10 distributed 
as part of HEASOFT (version 5.3). For each of the ob- 
servations, data were filtered using standard procedure 
given m the RXTE Cook Book@. For extraction of back- 
ground, model appropriate for bright source^! was used 
for the 2001 data since the source was in high state dur- 
ing that period. For 2003 data, when the source was in 
low state, the model appropriate for faint source^ was 
used. The resulting X-ray light curves for the two states 
of Mkn421 are shown in figure [4j 

Contemporaneous optical/NIR and radio data are 
also available for Mkn421 corresponding to the two data 
sets mentioned above. Optical data in V-band for year 
2001 was taken by WEBT (Whole Earth Blazar Tele- 



http://heasarc.gsfc.nasa.gov/W3Browse 



see http: / /heasarc. gsfc.nasa.gov / docs /software /lheasoft 
http://heasarc.gsfc.nasa.gov/docs/xte/recipes/cook_book.html 

(pca_bkgd_cmbrightvle_eMv20031123.mdl) 

(pca_bkgd_cmfaintl7_eMv20030330.mdl) 
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scope) using KVA-telescope on La Palma [7 . For 2003, 
Near Infra Red data in J-band was taken at Gurushikhar 
observatory, Mount Abu [8]. Radio data was taken by 
Metsahovi radio telescope at 22 GHz and 37 GHz [5]. 



4 Results 

Most of our observations pertain to Mkn 421. The upper 
panel of figure [5] shows the daily average of TeV 7-ray 
rate as obtained from PACT for Mkn 421. In the lower 
panel of this figure the daily average of X-ray photon rate 
as obtained by ASM (All Sky Monitor) on board RXTE 
is shown for comparison. PACT sensitivity is such that 
to have resonable signal to noise ratio we need to add 
data from several nights. Therefore we have obtained the 
time averaged flux of TeV 7-rays from Mkn421 by com- 
bining all data obtained using PACT during 2000-2005. 
This is found to be 4.5±1.9 xl0~ 12 photons cm~ 2 s _1 
above 1.2 TeV. Since Mkn 421 is observed at an average 
angle of about 30° with respect to zenith, energy thresh- 
old of PACT for these observations corresponds to 1.2 
TeV. In terms of Crab units this flux is 0.3 units. The 
overall significance is low due to systematic errors and 
is 2.3 a. This integral flux from Mkn 421 as obtained by 
PACT is shown in figure [51 The flux of 7-rays during the 
most intense flare and lowest activity state as obtained 
by Whipple group is also shown in this figure. The dot- 
ted line represents the quiescent flux measured in 1995 
[TO] and solid line represents the flux measured during a 
flaring state in 2001 [TTJ . During the period of our obser- 
vations, there were a few flares in 2000, 2001 and 2004. 
So our time averaged flux is expected to be between the 
two extreme limits given by Whipple and other Groups. 

We have not seen any significant 7-ray flux from other 
blazars we have observed (Mkn501, 1ES1426+428 and 
ON231). Earlier, in 1997, a huge flare was detected from 
Mkn501 12. 13lfl4]. But during our observation period 
this source was in low state. It was much weaker than 
Crab Nebula. We have estimated 3cr upper limit on 7-ray 
flux from Mkn 501, as 1.22X10 -11 photons cm" 2 s" 1 (> 
1.2 TeV). Average zenith angle for these observations was 
about 30°. This upper limit corresponds to 0.75 Crabs. 
This upper limit is shown in figure[71 The dotted line rep- 
resents the quiescent flux measured in 1995 [JO] and solid 
line represents the flux measured during the flaring state 
in 1997 by the Whipple group [T5]. Blazar 1ES1426+428 
is a distant object at z of 0.129. Whipple and HEGRA 
groups have detected TeV 7-rays from this source with 
long duration observations [T5I|16| . PACT is less sensi- 
tive than Whipple telescope. Also, we do not have very 
long coverage for this source. We estimate 3<r upper limit 
on 7-ray flux for 1ES1426+428 to be 1.34X10" 11 pho- 
tons cm~ 2 s _1 (> 1.2 TeV). Here energy threshold is 
in accordance with average zenith angle of 30° during 
these observations and upper limit corresponds to about 
0.8 Crab flux. ON231 is a LBL type blazar and till date 
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Fig. 5 Average 7-ray rate per night from 2000 to 2005 from 
PACT (upper panel) along with daily average X-ray photon 
rate obtained by ASM (All Sky Monitor) on board RXTE 
for Mkn 421. PACT data points include statistical as well as 
systematic errors. 
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Fig. 6 Integral energy spectrum of Mkn421. The data 
point (asterisk) with error bar represents the time averaged 
integral flux of Mkn 421 obtained from PACT observations 
during 2000 to 2005. Dashed line represents the flux mea- 
sured by Whipple during a low activity state and solid line 
represents their flux measurement during 2001 flaring state. 



TeV 7-rays are detected only from HBL blazars. But, ob- 
servations by EGRET on board Compton Gamma Ray 
Observatory (CGRO) have shown a hard power law en- 
ergy spectrum (photon spectral index a = 1.73 ± 0.18) 
extending upto about 10 GeV with no sign of any cut- 
off. Because of its hard spectrum it was thought that 
ON231 may be detected at higher energies and hence is 
a potential TeV 7-ray source. GeV/TeV 7-rays have not 
yet been detected from ON231 so far jl7U18j . ON231 was 
observed by PACT at an average zenith angle of about 
10° with energy threshold about 800 GeV. We have es- 
timated 3tr upper limit on 7-ray flux from this source 
as 2.50x10"" photons cm" 2 s" 1 (> 800 GeV) which 
corresponds to 0.83 Crab units and is shown in figure [5] 
along with other results. 
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Mkn 501 Spectrum 




Fig. 7 Integral energy spectrum of Mkn 501. The down 
arrow represents the 3<r upper limit on the flux of TeV 7-rays 
from this source obtained by PACT. Dotted line represents 
the flux measured by Whipple during a low activity state 
and solid line represents their flux measurement during 1997 
flaring state. 
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Fig. 8 Upper limits on the flux of GeV/TeV 7-rays from 
ON231 as obtained by Whipple, STACEE and PACT groups. 
EGRET spectrum is shown by solid line, which is extended 
upto TeV energies using dotted line. 



4.1 Spectral Energy Distribution of Mkn 421 

We have derived the spectral energy distribution of Mkn 
421 during high and low states. The former is based on 
contemporaneous radio, optical, X-ray and 7-ray obser- 
vations while the later is based on contemporaneous ra- 
dio, NIR, X-ray and 7-ray data. 

Spectral analysis of X-ray data was done using XSPEC. 
Spectral data from both the data sets were fitted by cut- 
off power law with line of sight absorption. Line of sight 
absorption was fixed to neutral hydrogen column den- 
sity at 1.38 x 10 20 cm 2 %F5\. The best fit photon indices 
for 2001 and 2003 data are found to be 2.05±0.03 and 
2.40±0.03 and cutoff energies are about 24.9±0.26 keV 
and 23.9±2.4 ke V respectively. 

Figure M shows spectral energy distributions (SEDs) 
obtained for Mkn 421 involving multiwaveband data for 
2001 and 2003. We have used the time averaged flux of 
7-rays for Mkn 421 in the SED plot of 2003. During the 
multiwavelength campaign in 2003 this source was at 



quiescent state as mentioned earlier. Mkn421 was at low 
state during most of the PACT observations from 2000 
to 2005 and hence the time averaged flux is expected to 
be closer to the quiescent state flux of Mkn421. For 2001 
flare we have selected PACT data overlapping with X- 
rays. After applying selection cuts we are left with 4.8 
hours of data which is simultaneous with X-ray data and 
the corresponding 7— ray flux is 3.4±1.4 x 10~ n photons 
cm~ 2 s _1 . We have fitted the SEDs of both these states 
with a simple one-zone SSC model (for detailed descrip- 
tion of the code see Krawzynski et al. [20]) as shown in 
figured solid line for 2001 and dotted line for 2003. This 
model assumes spherical blob of radius R and uniform 
magnetic field B, moving with respect to the observer 
with the Doppler Factor S, which is filled with a homo- 
geneous non-thermal electron population. Fit to SED of 
Mkn 421 spanning X-ray and 7-ray energies using SSC 
model has been attempted by number of authors for low 
as well as high state of the source. The fitted parameters 
broadly fall under two categories : some prefered larger 
Doppler factor 8 of about 50 ([21], [22] and [23]) while 
some others used lower value of S in the range 10-20 [24] , 
We have tried both these cases (5=50 and 14) and the 
fits for 6=50 are shown in figure [5J There is substan- 
tial change in magnetic field (B) and radius of blob (R) 
in the two states. We obtained B(magnetic field)=0.20 
G, R(radius)=2.3xl0 13 m and w(electron density)=0.03 
erg/cm 3 for the low state in 2003 with <5=50. For the high 
state in 2001 the respective values are 0.57 G, 1.3xl0 13 
m and 0.1 erg/cm 3 . For <5=14 the values for B, R and w 
are 0.28 G, lxlO 14 and 0.01 erg/cm 3 for the low state 
and 0.4 G, 9x 10 13 m and 0.04 erg/cm 3 for the high state. 
However it is also possible to fit the SEDs by varying R 
and B but keeping 6 same. 

Synchroton peak in the SED for 2001 is located at 
higher energies compared to that of 2003 as shown in 
figure [H implying spectrum hardens as flux increases. 
Even though one zone SSC model fits high energy emis- 
sion nicely (for 2001 and 2003) it fails to take into ac- 
count radio and optical fluxes. If we assume additional 
electron populations, as suggested by Krawzynski et al. 
for 1ES 1959+650, responsible for the low energy emis- 
sion, then the SED's of these components could also be 
fitted. An extensive study of Mkn421 was carried out by 
Blazejowski et al. [24] for the period 2003-2004 involv- 
ing radio, optical, X-ray and 7-ray data. There is good 
agreement between SSC parameters obtained by them 
and here using 5 ~ 14. 



5 High Altitude GAmma Ray observatory 

A 7 telescope array, called High Altitude GAmma Ray 
observatory (HAGAR), is being built at Hanle in the 
Himalayas, at an altitude of about 4.3 km, above mean 
sea level [25] . HAGAR is based on wavefront sampling 
technique like PACT. These 7 telescopes will be in the 
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SED of Mkn 421 




!/[Hz] 



Fig. 9 SED of Mkn 421 for 2001 (flare state) and 2003 (qui- 
escent state). SSC fits for both these states are shown by solid 
and dotted lines respectively. Three one-zone SSC models are 
used to fit X-ray, optical and radio data in each case. 




Fig. 10 Photograph of one of the telescope installed at Hanle 



form of a hexagon with an intcrtclcscope spacing of 50 
m. Each telescope has 7 para-axially mounted parabolic 
mirrors of diameter ^0.9 m with a photomultiplier tube 
at the focus of each mirror. The atmospheric attenua- 
tion of Cherenkov photons at Hanle altitude is ~ 14% 
as compared to ~ 50% at sea level. The Cherenkov pho- 
ton density near the shower core at Hanle is higher by 
a factor 4-5 compared to that at the sea level for show- 
ers of same energy. These features effectively reduce the 
energy threshold of HAGAR which is estimated to be ~ 
60 GeV for vertically incident 7-ray showers. Sensitivity 
of HAGAR would be such that it will detect the Crab 
at 5<7 level without any hadron rejection in ~ 2 hours. 
Figure [10] shows the photograph of the first telescope 
commissioned in June, 2005. Comissioning of remaining 
telescopes is underway. All 7 telescopes are expected to 
be operational by middle of 2007. 



6 Conclusions 

Out of 4 blazars which we have observed using PACT 
only Mkn 421 was reported (by other experiments) to 
be in flaring state on few occassions between 2000-2005. 
We have estimated average integral flux for this source 
by combining all the data from 2000 to 2005. For others, 
Mkn 501, 1ES1426+428 and ON231 we have given 3<r 
upper limit on 7-ray flux. 

We have found from the study of SEDs of Mkn421 
that X-ray and 7-ray emissions during high and low 
states are correlated. Synchrotron peak of SED of 2001 
is at higher energy compared to that of 2003, suggesting 
that peak shifts towards higher energy as flux increases. 
There are significant changes in SSC parameters for these 
two data sets at X-ray and 7-ray energies. One zone SSC 
model can not fit all data, introduction of additional zone 
improves the fit at lower energies. But there is almost no 
change for SSC parameters at lower energies during qui- 
escent and flaring state of Mkn421. 
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